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£\j ' Abstract 



We give a new method for the reduction of tensor integrals to finite integral repre- 
sentations and UV divergent analytic expressions. This includes a new method for the 
handling of the 7-algebra. 



In calculating two loop corrections to the Standard Model one is confronted with two main 
problems. One is the analytical difficulty of integrals involving different masses. Often one 
is restricted to approximations as for example zero mass assumptions for light particles |Q . 
C"| , A further problem stems from the difficulty in the algebraic sector, where the increasing 

number of terms for the tensor structure has to be handled. For the case of two-point 
functions an elegant method is available to reduce all tensor integrals to a basic set of scalar 
■ integrals ||] . Nevertheless it is clear that these methods will run into difficulties when applied 

to three- or four-point functions. Also, for arbitrary mass cases, one is restricted to either 
numerical evaluations or asymptotic expansions || . This is due to the fact that for arbitrary 
mass cases it is not possible to express the generic scalar integrals through known special 
functions pj. 

In this paper we will present a general method to express two-loop integrals in terms of 
finite integrals suitable to numerical evaluations plus a set of products of one-loop integrals 
containing the UV singular part. The method presented here is as well applicable to arbitrary 
tensor integrals but will prove especially powerful when applied to two-loop graphs directly, 
as it will bypass the whole tensor structure problem by constructing what can be called the 
characteristic polynomial of the graph. 

The final result is appropriate for the use of integral representation generated from the 
ones in || but the method can also be applied with other choices for the generic integrals, 
as the way how we handle the tensor structure is independent from the analytical approach 
one uses. 

The methods presented in the following apply to two-loop n-point Green's function for 
arbitrary n, but in examples we will restrict ourself to two-loop two- and three-point prob- 
lems. We will demonstrate the use of our method on the topology of the following figure. 
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Fig.l: The master two-loop two-point function; the labels denote the prop- 



A general two-loop integral has the following representation 

J>(ni,nfe) 

''NiMNk' 



j d n l,l"k- ■ (I) 



where (ni, n k ) denotes the rank of its tensor structure in I, k loop momenta. Here Ni is the 
part of the denominator containing only propagators depending on I, Nk the same for loop 
momentum k and M contains the propagators involving both loop momenta. 

We omit in the following the cases where terms in the numerator cancel the M part of 
the denominator completely out as these terms give rise to products of one-loop integrals 
only. 

Let us assume I has an overall degree of divergence ui and degrees of divergence loi , u k 
for the subdivergences. To make this integral convergent we have to subtract its I resp. k 
subdivergent behaviour as well as its overall divergence. 

To this end let us define the following 

Ni = Ni | TO . =0j9i= o, 
Nk — Nk | mi =o,9i=o, 

M — M | mi=0 , 9i =o, (2) 
(NiM - NiM) jl 



V 1 = 



{NiM)i> 
(N k M - N k My* 



where \ mi =o,qi=o means to evaluate the corresponding expression with all masses and external 
momenta set to zero. 

Now consider the replacement 

m Ni ' 

which improves w and uii at least by ji and leaves Uk unchanged, and a similar replacement 
for k. Consequently 

1 K j *L>\ (3) 



NiMNk NiMNk 
results (at least) in an improvement 

u — » u> + ji + jk = 
u>l —> LOI+ ji = 

LU k —> LOk + jk = 



u'k- 



Let us introduce the notation oS > {co' > 0,u[ > 0,uj' k > 0}. Let us further define the 
notation Co := C — 1 for arbitrary expressions C so that 



Til _ . 

L o - 



(NiM)oi 



31 



1=1 



J2( J t)(- N i M y(tiiMy- 
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K 3k 



1 



[L jl K jk ] 



{N k M)i* 
1 



1 



(NiM)i' (N k M)i k 
3i\ fjk 



]T ( J *)(-N k M) k (N k My*- 1 

ji jk 

EE 



l \k 



1 = k = 

l+k>0 



{-NiM)\NiM)^- l {-N k M) K {N k M) 



(4) 
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where we emphasize that the sums start with I = 1 resp. k = 1. 
We then have the algebraic identity 

rp(ni,n k ) T^(«i,n fc ) rp( niy n k ) 

= \L jl K jh ] \L j 'K jk ] n Vji,j k . 

NiMN k NiMN k 1 J NiMN k 1 Jo J ,J 

Now there clearly exists for a given T and given Ni , M, N k powers ji , j k such that 
u>' , u>i, u)' k are all positive, u' > 0, so that the first expression on the right hand side of 
equation (^) is finite. So it can be calculated in D = 4 dimensions. But by inspecting the 

f2( J !)(-NiM) l (NiMy- 1 , 

E ft) (-N k M) k (N k My»- 



sums 



L Jl = 



1 
1 



(JV fc M)i* 



Lfc=0 



we can identify all terms in the expansion which fulfil lo > 0. As these terms can be 
calculated in four dimensions they cancel against the corresponding terms of the second 
expression in the right hand side of equation (||). Note that to identify the divergent part a 
simple counting of dimensions of masses and exterior momenta in the numerator of equation 
(||) is sufficient. 

So we define the following expressions 



tJi 
L 



^0 



(NiM)i 



(N k M)i» 



(NiM)i 



1 

(N k M) jk 



1 = 

w>0 



h -i 



k=0 

w>0 



( Jk h )(-N k M) k {N k My k - 



1 = 1 
w>0 



3k 



£ ( 3 l)(-NiM) l (NiM) 



k = l 

w>0 



J2 ( J !;)(-N k M) k (N k My k 



-i 



(6) 



where X)w>o means summation only over divergent contributions. We obtain 
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NiMN k 
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N t MN k 
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Still the first of the two terms on the right hand side of equation (0) can be calculated 
in D = 4 dimensions. We will now consider an example which will be useful to demonstrate 
the method. Consider Fig.(l) with propagators 



Pi = I 2 - to 2 



l! 

P 2 = (l + qf-ml, 
P 3 = (l + k) 2 -m 2 3 , 



P 4 = (k-q) 2 -m 2 , 
P 5 = k 2 -m 2 5 , 



and choose 



y(2,l) 


— Ilk 


Ni 


= PlP-2, 


M 


= Ps, 


N k 


= PaP 5 , 


N, 


= l\ 


M 


= Q + kf, 


N k 


= k\ 


u = -l, 


u>l = 0, u k = 1, ji 



so that we have 



2, Jfc=0, 



T(2,l) T(2,l) T(2,l) T^'^NiM 

-L 2 = 2-^- 



NiMN k NiMN k NiMN k (NiM) 2 N k 

-Lr, — —2— =- 



NiMN k NiMN k {NiM) 2 N k 

^(2,1) . y(2,i) r (2,i) T { - 2 - 1 \2l-q) 

-L 2 = 



NiMN k NMN k ^MNk N[Ml 2 N k 

T (24) /2 T (2,i) | T( 2 ' x )(2/-g) 



£2 = ____ + . . v (8) 

iV ; MiV fc NiMN k NiMl 2 N k 

where we used that Pi = Pi — to 2 , P2 = P2 + 21 ■ q + q 2 — rn 2 , P3 = P3 — m§. Pi means again 
the suppression of all masses and exterior momenta in the corresponding propagator. 

We will continue later with this example and discuss now how to calculate the UV 
divergent part of equation (Q) , which is contained in the second term of the right hand side 
of this equation. These integrals are of the form 



rp(m,n k ) 



{i 2 y{i + k)Wk 



(9) 
(10) 



(k 2 y(i + kyNi' 

0. (11) 



(i 2 Y(i + ky(k 2 ) r 



The second equation, (|l0|), is just a relabeling of (||). 

The third equation, (p_l[), vanishes identically in dimensional regularization. The integer 
powers i,j in the above equations are determined by ji,jk- 

Let us comment on the vanishing of the integral in ( |TT| ) here. It might happen that this 
cancellation is due to a mutual cancellation of an infrared and an ultraviolet divergence, as it 
appears in the integral J d D l /I 4 . Now the UV singular part of the integral in (|iT| ) will have 
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its counterpart in the final integral representation in (Q), where it is needed to cancel some 
other UV singular parts. So it seems that we have installed a spurious infrared singularity in 
our integral representation. This should then be regularized by some appropriate IR cutoff 
and this would install a corresponding scale in the integral in ([ll]) which would make it non- 
vanishing. But these infrared singularities appear in the integral representation derived in 
H always as endpoint singularities. As this problem with spurious IR singularities can only 
appear when simultaneously n\ > and tih > it will (in the Feynman gauge) be restricted 
to the two-loop boson self-energies. Nevertheless it is true in general that IR singularities 
appear as endpoint singularities in the method of [j| which avoids parametrizations. Let us 
study these endpoint singularities a little bit more. Investigating the case 



Ni 


= Pi, 


M 


= P 3 , 


N k 


= Ps, 



we find the following expression for the subtraction of UV divergences 

1 1 1 M~m\ 



N t MN k N[MN k NiMN k N t M 2 N k 



(12) 



The last term ~ is the problematic one. The integral representation for it can be written 
as 



oo 



X 



' dy \x\\x\ + \y\ + \x-y\' ^ 



where we used yj x 2 + irj — \x \ in the limit r/ — > 0. Equation ( |l3| ) has an apparent endpoint 
singularity at x — 0. But these endpoint singularities are a well defined distribution in the 
ry-limit when interpreted as Hadamard's Finite Part HF P ||. It s relation to the propagator 
is known for long 0. In fact, the distribution dxj 'yj x 2 + ir\ equals HFP[1/ \x\]. Using 
this, we handle these endpoint singularities with the help of HFP and find the corresponding 
integral representation 

f°° 1 1 
dx / dy(-, r-j — ; j — ; ; r + X <-> — X) + . . . . 



Here the dots correspond to finite terms for which the endpoint singularity in x has canceled 
out after the y integration and so HFP reduces to an ordinary Riemann integral. As HFP 
is a local operation which modifies test functions only at the location of the singularity our 
UV behaviour remains unchanged and the sum in equation (^) is still UV finite. This can 
also be explicitly checked by working out all the different cases for the modulus function in 
equation (|l^) in detail. So the integral representation for equation (|l2|) is free of IR and UV 
singularities. One can also use HFP to separate nonspurious IR singularities when they 
appear. As for the integral ([ll]) to appear we need UV singular behaviour in both loop 
variables, it is clear that the above subtlety can only arise when there is a logarithmic UV 
divergence in the second loop integration, which is exactly the case for the above example. 
In this sense the above example is generic and serves as a proof in general. 

We can still maintain the vanishing of equation ( |i"T| ) by using HFP in our integral 
representations. This can also be interpreted as a consistency between HFP and DR. 

Let us continue now our general discussion. We can do one integration in equation (|^) 
leading to remaining integrals of the form 

C [ d D k—^- , (14) 

where a is an integer only when D = 4 and T is some tensor in k. The coefficient C will 
contain a pole in (D — 4) in general. This and possible UV divergences of the k integration 
itself forbid to calculate these integrals in four dimensions. 



Calculating the above integrals fll4| ) for arbitrary D and a is possible but a difficult 
analytical task. But there is an algebraic approach to these integrals. We will first transform 
the UV divergent integrals in equation (fjj) via a partial integration in such a way that the 
first I integration is UV finite and its UV divergences are shifted to the k integration. Then 
we use a method similar to the method above to split the final k integration once more into a 
finite integration, analytically computable in an easy manner and UV divergent integrals of 
massive tadpole structure. As a consequence everything is expressible by algebraic methods 
in one-loop functions. 

Defining 

l(l ' n ' m) ■= I dDl nJ^'," 1 ^ . w e have 
J {l 2 ) n {{l + k) 2 ) m ' 

n _ f iDi d y M1 ■■■i IH , s 

= (D + i-2n- m)I^ n ' m ^ - mi^-i.m+i) + mjfe 2 J (i,n,m+i) j 

j(<,0,m) =7 (i,m,0) _ y-^ ( 16 ) 

As D + i > 2(n + m) for UV divergent integrals, we have (D + i — 2n — m) 0. Therefore 
we can always solve equation ( |f6| ) for /(^™^ m ) in the following. 

We can repeat this procedure until we have expressed j( l >"' m ) through integrals of the 
form ' m ' with powers m',n' such that they are all finite in the / integration. 

We will end up with integrals of the form 

j^fe^ii ■ ■ ■ k fJ/r (k 2 )^ 



d u k 



where again (3 is an integer only if D = 4 but all coefficients of these integrals are now finite 
in four dimensions. 

But the k integration is still UV singular, and to find an expression for these integrals 
we do a last replacement 

± = where (17) 

JO := and 

N k := N k |, 4=0) 

where it is understood that when expanding the numerators the sums run only over the 
divergent parts, as indicated by the " on the various K . Note that we have chosen propa- 
gators with vanishing exterior momenta but non- vanishing masses to profit from symmetric 
integration properties and we remind the reader that we cannot set the masses to zero be- 
cause this would result in vanishing tadpole integrals. So we end up with finite integrals 
and separated UV divergent integrals of the form 

k k (k 2 )P 

which are, as they are massive tadpole integrals, easy to evaluate in D dimensions. The 
finite integrals in ( |l7| ) are related to standard massive one-loop integrals with integer powers 
of propagators. 

Let us continue our example now. By inspecting equation (||) we see that we have to 
consider l( 2 > 2 <V and J^ 3 - 1 ) . The last one is already finite in the / integration. For the first 
one we have 



7(2,2,1) : ~ :l 



(D-3) 



Z fc 2r(2,l,3) _ u 2 r(2,2,2) 

(D - 2) 



(> 



The remaining integrals are of the form 

h h 



Dt, "'Ml ■ ' ' "7*< 



p 4 p 5 (k 2 y+— 

By using P 5 = P 5 ,P 4 = P 4 — 2k ■ q + q 2 we separate this according to the above equations 
dl7| ) . Note that as a consequence of our subtraction operation only terms with an even power 
of k in the numerator give contributions to UV divergent integrals. The finite integrals are 
of the form 

Dl k fil ...k tti (-2k-q + q 2 Y 



d u k 



p 4 p i p 5 (k 2 y+ 1 



where only terms of even power in k contribute in the numerator. 

Via partial fraction decompositions in k 2 ,P5,P 4 this can be reduced to standard one- 
loop integrals with integer powers of propagators which are algebraically related to one-loop 
integrals with unit powers of propagators via mass derivatives. 

All these reductions and separations can be handled by symbolic calculation routines. 
So it is possible to reduce an arbitrary two-loop tensor integral into finite integral repre- 
sentations and a standard set of one-loop integrations. For the examples discussed in this 
paper this was always possible using REDUCE § on a 486 notebook within less than 5 min- 
utes. Note that even for the two-loop four-point function Nh is not worse than a one-loop 
three-point function (apart from the trivial topology arising from self-energy insertions), so 
that the resulting one-loop integrals are relatively simple. The knowledge of arbitrary tensor 
integrals of one-loop two- and three-point functions with integer powers of propagators is 
indeed sufficient for the calculation of two-loop problems by this method. 

Having demonstrated the handling of the UV divergent part let us now discuss the 
integral representations for the finite integrals. In Q an integral representation is given 
which transforms the integral to a twofold integral over the parallel space variables. This 
can be generalized to arbitrary tensor integrals which appear as polynomials in parallel and 
orthogonal space variables in the numerator. So we have to integrate out the orthogonal 
space variables. This can be easily handled by use of partial integrations and applying the 
residue theorem afterwards. But a little modification arises for the terms which subtract the 
UV divergence. Here we cannot use a partial fraction decomposition in the Pi, Pi propagator 
as these propagators are equal. Doing the z and fcj_-integration as in j5) we are left with 



J dl± - (p r]2 y [ lo S terms ]; i = 2,3. 



But we can integrate these terms most easily by rewriting 

1 Id 1 1 



(i 2 - i\Y i\ dp {il m) 



(18) 



A similar remark applies to the three-point two-loop functions. The case of four-point 
functions will be discussed elsewhere, but will give no modifications in general ||. Never- 
theless we include the principal reduction of two-loop box integrals in the discussion below. 

Up to now we have studied particular tensor integrals only. Let us now discuss a method 
how to calculate tensor integrals more directly. We again use the notation of || for the 
parallel and orthogonal space variables and will discuss as a specific example the following 
graph, which appears in the flavour changing fermion self-energy. The two-loop radiative 
corrections to order 0(a s ) x g 2 are currently being investigated for these self-energies [fL0| . 



cles. 



-•4 q i Fig. 2: Again the master topology; the labels specify propagators and parti- 
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Let us decompose the Clifford algebra 

/ 

{7||j7-l} 

where we simply use the fact that we can split a loop momentum in DR in its component in 
the direction of the exterior momentum q = q\\ and its orthogonal complement. The same 
splitting can be done for higher dimensional parallel spaces (the linear span of the exterior 
momenta) . This gives rise to a corresponding splitting in the Clifford algebra which can also 
be written down in a covariant manner 

= 

l±n = 

7-Lm = 

7||p = 
9^v = 

{7|k>7-L/i} = 

The generalization to higher dimensional parallel spaces is clear. 

We can now use this decomposition for a reordering where we bring all parallel space 7 
matrices and the anticommuting 75 p~T| ] to the left so that we have 

for an i dimensional parallel space and a rank k tensor in the orthogonal space. The rii are 
the number of appearances of 75 and parallel space 7-matrices. Using 7s 2 = l;7n. = 1 we 
can simplify the above expession. For the case that we do not have free Lorentz indices in 
our Green function we can also simplify the orthogonal space part of this expression. In 
this case we are allowed to replace the orthogonal space j± matrices identically by metrical 
tensors g±, (gj_n = (D — i)), even when we do not apply a trace to the string. Note that 
this means that all the lengthy trace and 7 algebra calculations of perturbation theory 
will not appear for the case of no free Lorentz indices. It also means that we only have 
to consider contributions with an even number of 7-matrices in the orthogonal space. To 
prove this statement note that as long as we have no free indices (after contracting out all 
double indices) 7^-matrices can only be generated by the loop momenta /, Iji. As the only 
covariant available for the orthogonal space integration is the total symmetric combination 
of metric tensors g± (the Levi Civita tensor does not appear in symmetric integration) we 
only need the symmetric part of the Clifford algebra in orthogonal space. This corresponds 
to a replacement of 7j_ matrices by g± tensors, as all indices have to be contracted by a 
total symmetric tensor. 

This method can be generalized. In the case that we cannot avoid free indices we still 
can use the above argument for the 7^ matrices generated by orthogonal space components 
of loop momenta. But note that in most cases it is possible to contract free indices with 
metrical tensors or exterior momenta without losing information. For example, to determine 
the form factors appearing in radiative corrections of the fcrmionic vertex, one can contract 
with exterior momenta p,p' say. Then one can use all advantages of the case of no free 
indices. 



= l\\T-l±i7x, 
= fc^ 1 - fc_u7l, 
= 0, 



l-q 

n 2 



l-q 



aw 



9 q^,^fj, 9 q^ 

q z q A 



= 0. 
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Our example was calculated using this method by installing the proposed ordering algo- 
rithm via simple LET rules in REDUCE. Comparing this with the same calculation using 
the high energy packet in REDUCE, it was found that the method we suggest here was 
orders of magnitude quicker. Even for the one-loop case we can report on a calculation of 
two- and three-point functions in the unitary gauge in the SM where significant advantages 
were obtained for the 7 algebra Jl2^ ]. 

This method fails for the anomaly. The interesting covariant is ~ ^ il vpaV P \P%. an( i there 
exists no tensor constructed from symmetric metrical tensors and exterior momenta p\ , P2 
to give a nonvanishing contraction with this covariant. One nevertheless can calculate the 
resulting trace by this method. But then one has to consider the case of 75 with six 7 
matrices, two in the linear span of pi,P2, two in the orthogonal complement but still in 
four dimensional Minkowski space and two with arbitrary indices, which gives rise to the 
non-cyclic trace of Q. 

Continuing with our example we are left with the tensors 1,75,1^,75^. Each of these 
tensors is multiplied by an integral representation for the finite part, which is of the form 
(using (0),([l8|) and the results in §) 



1 



(Pi-P 2 )0P4-P E 

— [(c(3, 0) + c(2, 1) + c(2, 0) + c(l, l))£ sub ] U=0 - I ■ q-^ [(c(3, 0) + c(2, l))£» b ]|^ 
■ — [(c(0, 3) + c(l, 2) + c(0, 2))£ reb ]| M=G + k ■ q-^ [(c(0, 3) + c(l, 2))/C Bllb ] u=0 . 



[c25-^325 — Cl5^315 + Cl4-^314 _ C24.C324] 



We list the explicit coefficients C25, C15, C24, C14, c(i,j) in an appendix. For the case of Fig. (2) 
it took only about 20 seconds in REDUCE to generate them for the complete graph. 

Let us summarize the situation for the most interesting two-loop graphs. We assume 
that the corresponding calculations are done in a renormalizable gauge, though it is not 
necessary to restrict to this case. 

Two-point functions: 

We have two topologies. For the trivial self-energy insertions we can always arrange that 
power counting improves by 2n; (assuming I is the subloop momentum), as we can always 
arrange exterior momenta not to flow through the subloop. For this topology we then have 
rifc = always, while n; = 1 or ni = 2 corresponding to to — 1 or lu — 2, that is corresponding 
to boson or fermion self-energies. 

For the other (master) topology of our example we have (rife = 0, n; = 2) or vice versa if 
lu = 1 and (ni = 2, = 1) or vice versa if u> = 2. 

Three-point functions: 

We have three topologies, the ladder topology, the crossed topology and the self-energy in- 
sertion topology. For the last one we always have (n; = 1 , nk = 0) while the other ones 
reduce to (rik = 0, ri; = 1) or (rife = 0, ni = 2), depending what kind of vertex one considers 
(the fermionic one has only logarithmic overall divergence, while e.g. triple boson couplings 
are linear divergent). Note that for the crossed topology, as M contains two propagators, the 
UV divergent one-loop integrals reduce to two-point functions and so this more complicated 
topology is simpler in its divergent structure. This corresponds to the fact that it has no 
proper subdivergence, while the ladder topology has a vertex type subdivergence and thus 
involves three-point one-loop functions. 



Four-point functions: 

We need four exterior bosons to have a UV divergence of logarithmic degree. We then have 
u) = and the case (ni — l,nk — 0) is always sufficient. Again the complicated topologies 
involve simpler one- loop integrals. 



9 



We see that only for boson self-energies we have the case (n; / 0, rik ^ 0) so that only 
for them do we have to use HFP explicitly. Note also that all (m = 1,71* = 0) cases reduce 
to a very simple subtraction which simply sets the exterior momenta and masses to zero 
in the subtracted term, and this is sufficient especially for the two-loop fermionic vertex 
corrections. 

We conclude that an easy separation of divergences for two-loop functions can be achieved. 
We saw that for the subtraction of the UV divergences of the real two-loop part very simple 
one-loop functions were sufficient. Especially only terms of order 1/(D — 4) were used. This 
was not unexpected as the order l/(D — 4) 2 poles of a two-loop graph can be obtained from 
products of one-loop functions. This is true in general for the leading divergence in a n-loop 
calculation. So we will find the leading divergences in the terms where the M part of the 
denominator is cancelled, giving rise to decoupled one-loop integrals. For our example in 
Fig. (2) we list the corresponding expressions in the appendix. As the individual terms of 
equation (0) correspond to terms subtracting Z divergences, k divergences and overall diver- 
gences, one might wonder if it is possible to implement the whole renormalization program 
(that is the addition of one-loop counterterm graphs) at this level to obtain even simpler 
expressions. We intend to investigate this question in the future. 
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Appendix 



In this appendix we list the relevant polynomials in Iq, k^ for the final integral represen- 
tation of the graph of Fig. (2) explicitly. Our notation follows ||. We list these results to 
give an idea of the complexity resp. simplicity of the integral representation for a complete 
graph. We see that the integral representation, compared with the pure integral represen- 
tation for the scalar master function, is similar in its qualitative and quantitative difficulty. 
Let us first list the propagators of Fig. (2): 

P 1 :=l 2 -m 2 j P 2 :=(l-q) 2 , 
P 3 := (l + k) 2 -m 2 , 
Pi := (k + qf - mj P 5 := k 2 - m 2 . 

We define q := \[q* and := <j[/q. 
ci5 ■= 8 • 75^ • (- (2 • q ■ l 2 ) 

- 2 • q ■ l ■ k + m 2 • Z + m 2 w ■ k - m 2 • k + m 2 ■ l + m 2 ■ k ) 
+ 4 • 75 • n%j ■ (- (2 • q ■ Z ) - 2 • q ■ k - m 2 w + 3 • m 2 + to 2 ) 

+ 8 • 4 ■ (- (2 • q ■ Z 2 ) - 2 • q ■ Z • k a + m 2 w ■ Z + m 2 w ■ k - m 2 ■ k 

+ m 2 • Z + m 2 ■ Zc ) + 4 • mj ■ (2 ■ q ■ Z + 2 • q ■ k + m w — 3 • to 2 - to 2 ) 

ci4 := 8 • 75^ ■ (- (q 2 ■ l ) 

— q 2 ■ k — 2 • q ■ Zq — 4 • q ■ Z • k - 2 • q ■ k 2 + m i ■ Z + to 2 • l 
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+ m 2 ■ ko) + 4 • 75 • rrij ■ (q 2 - 2 • q ■ l + 2 • mf + mf) 

+ 8 • 4 • (- (q 2 ■ Iq) - q 2 ■ k - 2 • q ■ I 2 - 4 • q ■ l ■ k - 2 • q ■ k 2 

+ mf ■ l + m 2 ■ l + m 2 ■ k ) + 4 • rrij ■ (-q 2 + 2 • q ■ l — 2 • mf — mf) 

c 2 5 := 8 • 75^ • (- (q 2 ■ lo) - q 2 ■ k 

+ mf ■ l + m 2 w ■ k - mf ■ k ) + 4 • 75 • rrij ■ (-q 2 - 2-q-k a -m 2 w 

+ 3 • mf) + 8 • }■ (- (q 2 ■ l ) - q 2 ■ k + m 2 w ■ l 

+ m 2 w ■ k - mf ■ k ) + 4 • rrij ■ (q 2 + 2 • q ■ k + m 2 w - 3 • mf) 

c 24 := 8 • 75^ • (- (2 • q 2 ■ l ) 

— 2 • q 2 ■ ko - 2 • q ■ l ■ k - 2 • q ■ fcg + mf ■ l ) + 8 • mf ■ rrij ■ 75 

+ 8 • i • (- (2 • q 2 ■ lo) - 2 • q 2 ■ ko - 2 • q ■ l ■ k - 2 • q ■ k 2 + mf ■ l ) - 8 • mf ■ m, 

c(3,0) :=8-/ - M- 75^ + 8 ■ l ■ M ' 4 
c(2, 1) := 8 • fc • 11 ■ j 5 4 + 8 • k Q ■ \i ■ 4 

c(l, 2) := 8 • lo ■ n ■ ibi + 8 • l ■ A* ' 4 
c(0, 3) := 8 • ko ■ 11 ■ 75^ + 8 • fc • \i ■ 4 
c(2, 0) := - (l6 ■ q-l 2 - 75^ + 4 • rrij ■ fj, • 75 - 16 • q ■ 1% ■ \ - 4 • mj ■ 11 

c(l, 1) := - (l6 • q ■ l ■ k • 75^) - 16 • q ■ l ■ k ■ \ 
c(0, 2) := — (4 • rrij ■ 11 ■ 75) + 4 • rrij ■ \i 



- In y^i+^-ml + k 2 + I Jo + fco|)) /{Pi - P 5 ) 

/C sub := (- In (^Jk 2 -ri + yj ' -m 2 + l 2 + \l + k 

+ In f^k 2 - I lo - q I + I ko + Jo |)) /(A - P 2 ) 

<X> := - (8 • 75^ • {lo + fco)) - 4 • mj • 75 - 8 ■ ^ ■ (Z + fco) + 4 ■ m^- 

where we used the symbol CO to denote the term where the M part of the denominator is 
cancelled, so that is corresponds to a squared one-loop graph. 
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